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We propose a general strategy to determine universal properties induced by a nearby phase tran- 
sition on a non-order parameter field. We use a general renormalizable Lagrangian, which contains 
the order parameter and a non-order parameter field, and respects all the symmetries present. We 
investigate the case in which the order parameter field depends only on space coordinates, and the 
case in which this field is also time dependent. We find that the spatial correlator of the non-order 
parameter field, in both cases, is infrared dominated and can be used to determine properties of 
the phase transition. We predict a universal behavior for the screening mass of a generic singlet 
■ field, and show how to extract relevant information from such a quantity. We also demonstrate that 

' the pole mass of the non-order parameter field is not infrared sensitive. Our results can be applied 

to any continuous phase transition. As an example, we consider the deconfining transition in pure 
Yang-Mills theory, and show that our findings are supported by lattice data. Our analysis suggests 
that monitoring the spatial correlators of different hadron species, more specifically the derivatives 
of these, provides an efficient and sufficient way to experimentally uncover the deconfining phase 
transition and its features. 
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I. INTRODUCTION 
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The relevant properties of phase transitions are best investigated using order parameters. However, sometimes it 
is profitable to determine new universal features associated with non-order parameter fields 0,0- Besides, in nature 
most fields are non-order parameter ones. To derive the general properties of a generic singlet field we couple this to 
the order parameter Q. In this way we can directly study the transfer of information, envisioned in 0, of the phase 
. transition properties from the order parameter to the singlet field. 

We analyze two cases: time dependent order parameter fields and time independent ones. The non-order parameter 
Oh, field, always time dependent, is a scalar under Lorentz transformations. Also, it is a singlet under any symmetry 
group, in particular under the symmetry group whose breaking is monitored by the order parameter. We further 
D \ assume the non-order parameter field to have a large mass with respect to all the other scales in play, and hence, an 
• associated small correlation length near the phase transition. Although time independent order parameter fields carry 
| the physical information about the phase transition, they do not propagate and cannot be canonically quantized. An 
example is the Polyakov loop, which is considered to be the order parameter of the pure Yang-Mills theory, and by 
construction is a function only of space. A generic singlet field in the Yang-Mills theory is the glueball, which is also 
a physical state of the theory. 

Time dependent order parameter fields are generally associated with physical states. These can be either composite, 
such as the chiral condensate, the order parameter related to chiral symmetry in QCD, or elementary, such as the 
Higgs field might be for the electroweak phase transition. Using general field theoretical arguments we demonstrate 
that the screening mass associated with the spatial two-point function of the singlet field is heavily affected by the 
nearby phase transition for both time dependent and time independent order parameters. The important fact is that 
we can predict the general behavior of the static two-point function of the singlet field, and associate it uniquely 
with the specific character of the phase transition. More specifically, the screening mass at one loop has a drop near 
the phase transition. This is due to the three dimensional nature of the screening mass, which makes it particularly 
sensitive to infrared physics, and it explains the transfer of information between the order parameter and the singlet 
field. When going beyond the one loop approximation the drop of the screening mass of the scalar singlet field, both 
for time independent and time dependent order parameters, is finite at the phase transition. Furthermore, the drop 
itself is controlled by the ratio of the square of the relevant coupling of the singlet field to the order parameter and the 
coupling governing the self-interaction of the order parameter. The result for the three dimensional order parameter 
has already been presented in 0. Here we provide details of the basic computations used in Q, especially for the 
higher order corrections, as well as a proof of the result for the screening mass in the time dependent case. For the 
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latter we also present the effects of the nearby phase transition on the pole mass of the singlet field. We demonstrate, 
somewhat surprisingly, that although we still observe a small drop, this physical quantity is not infrared dominated. 
So in general, only the spatial correlation lengths feel the presence of a nearby phase transition. 

The paper is organized as follows: In section we introduce the theory which is investigated in detail in section 
11111 for the case of an order parameter field that depends only on the space coordinates, and in section IIVI for an 
order parameter field that depends both on time and space. The information conveyed in section ITTT1 complements 
and enlarges the one presented in ref. Q. In section E] we conclude. 

Our main conclusion in both cases, space and space-time dependent order parameter fields, is that the information 
about the phase transition, encoded in the behavior of the order parameter field is transferred to, and obtainable 
from the singlet field(s) present in the theory. The present analysis suggests that monitoring the spatial correlators 
even of heavy hadrons provides an efficient and sufficient experimental way to uncover the existence and features of 
the chiral/deconfining phase transition. So if our results are phcnomcnologically applicable in relativistic heavy ion 
collisions at RHIC, they would lead to a clear signal for the existence of a deconfined phase at the early instant of 
such a collision. 



II. GENERAL SET UP 



We consider a temperature driven phase transition and work in a regime close to the phase transition. In order 
for our results to be as universal as possible, we use a renormalizable Lagrangian containing a field neutral under the 
global symmetries, and the order parameter field, charged under the symmetries, as well as their interactions. The 
protagonists of our theory are two real fields, h and x- The field h is a scalar singlet, while x transforms according to 
X — > z x with z S Zjy . 

While the generalization to Zjq is straightforward, we consider explicitly the case of ^2, which is suitable for 
understanding the deconfining phase transition of two color Yang-Mills. This has been heavily studied via lattice 
simulations [ijQ. The most general renormalizable potential is: 

V{h,x) = —h +9oh+ — h X +j^X + y h + 4\ h ' W 

The coefficients are all real and stability requires A > and t/4 > 0. We further assume <?i > and go < 0. At this level 
we did not yet commit on the space-time dependence of the fields, although we have normalized them as if they were 
living in 3 + 1 dimensions. This renormalizable potential can be considered, for example, as a truncation of the one 
presented in and can be used to determine some of the space-time independent properties of the vacuum. In order 
to go beyond constant field approximation, we must specify the kinetic terms. Since the singlet field h is physical, 
it has the ordinary standard four dimensional kinetic term, i.e. d^hd^h. For the x field we first consider the case in 
which x is time independent. In this case the associated kinetic term reads as VxV%. We will then consider the case in 
which x experiences also time dimension, and thus the kinetic term is d^x^X,- The time and space independent order 
parameter field, and its interactions with a singlet field were introduced in [lj , in order to understand the transfer of 
information from Polyakov loops to glueballs in pure Yang-Mills theories. The first studies of a renormalizable version 
of the previous idea, in which the order parameter is a function only of the space dimensions, were performed in 0. 
Here we review and extend the analysis of |2J, and show the fundamental differences and similarities with respect to 
a time dependent order parameter field. Before proceeding we state the assumptions under which we conduct our 
analysis for the time (in)dependent order parameter: i) The x field is light close to the transition, hence it dominates 
the dynamics, ii) The h field is heavy, and thus we can neglect its quantum and Boltzman suppressed thermal loop 
corrections. The extremum of the linearized potential (in the h field) is for 
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(X) 2 = -6^ with ml~ml x + gi {h) , and {h) ~ - ^(x) 2 ■ (2) 

Here m x vanishes at the phase transition. Having g\ > and go < 0, together with 3g 2 < Am 2 , assures the positivity 
of these expectation values, and with this choice also the extremum of the potential is a minimum. A more complete 
treatment would require to go beyond the linearized approximation. The temperature dependence of (x) and (h) is 
qualitatively sketched in figure ^ Note, that in Q the term linear in h is removed, at tree level, by shifting the field 
as follows: 



h-> (h) + h 



(3) 
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FIG. 1: Schematic behavior of the expectation values of the order parameter (x) (solid-line) and the non-order parameter (h) 
(dashed-line) fields close to the phase transition as a function of the temperature. Note that (h) is never zero and the constancy 
in the unbroken phase is due to the tree level approximation. 



Such term, however, is regenerated at higher orders, correcting the h expectation value. The effect of the shift leads 
to modified coefficients in the Lagrangian. For example: 

9a -> , 3i -> si + 32 (h) , ml x -> m\ x + g x (h) + y (h) 2 . (4) 

In the following we consider the fluctuations of h around its vacuum expectation value. 



III. TIME-INDEPENDENT ORDER PARAMETER FIELD : PROBING STATIC PROPERTIES 



The first difference between a purely spatial order parameter and a time dependent one is that in the first case no 
Bose-Einstein distribution emerges for x i n thermal equilibrium. More specifically, we are allowed by construction to 
consider only spatial fluctuations. Besides, since we assumed the field h to be heavy compared to the relevant scales 
at the phase transition, including the temperature, the h induced thermal corrections can be safely neglected. We 
postulate, as customary, that m x is a function of temperature and vanishing at the transition point. Interestingly, 
since the order parameter field is purely three dimensional, we can freely choose from which direction the symmetry is 
restored. Higher order temperature corrections cannot reverse the direction of symmetry restoration. As we discuss 
in detail in the next section, this is not the case if the order parameter field is time dependent. Having in mind 
the Yang-Mills deconfining phase transition, we assume that the Zi symmetry is restored at low temperatures. Our 
analysis is not affected when reversing the direction of the transition. 

We have chosen \ = x( x )j and decomposed the four dimensional field h into its Matsubara modes. After integrating 
over time the action reduces to an effective three dimensional one. In the three dimensional case the kinetic term 
for h and its self-interaction terms receive contributions from all Matsubara modes. However, only the zero mode 
contributes to the /ix 2 -interaction, which is the one driving the dynamics of the h field close to the phase transition 
0, as the x field becomes light with respect to h. Hence, we confine our discussion here to the theory which features 
the fields x and ho, and study the spatial fluctuations. For simplicity, in this section we denote ho by h which is taken 
to be directly the fluctuation field around its tree level vacuum expectation value. The three dimensional Lagrangian 
reads: 



C 3 = ±VhVh + iv X V X + \m 2 h 2 + \ml X 2 + T^f + VT 9 -± h X 2 + t| h\ 2 + Vf^h 3 + T 9 -±h* 



(5) 



where the coupling constants have the same mass dimension of the corresponding four dimensional theory. All of 
these are real with A > and 34 > . 

Using the renormalizable theory defined by the above Lagrangian, we now compute all one-loop corrections. We 
find the following contributions to the h one point function: 



o 



T 

— (gim x + g 3 m), (6) 

07T 
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which leads to temperature dependent corrections to the expectation value of (h). We turn our attention to the h 
two point function. The full expression at one-loop level is given by the following set of diagrams: 



= 2-0- + 2-Q- + e 9 

+ 12 — — + 18— Q— +18 



(7) 



Double lines indicate \ fields, while single lines stand for the h field. The number in front of each diagram is the 
associated combinatorial factor. We regularize the ultraviolet divergent contribution by subtracting the divergent 
terms. The second diagram is infrared divergent for T — ► T c and gives the dominant one loop contribution to the 
screening mass of h. In the limit of zero external momentum: 



-o- 



01 2 f d 3 k 1 _ g\ 

J (2tt) 3 (k 2 + m 2 f ~ T 32nm x ' (8) 



It is worth mentioning that (h(x)h(0)) is associated to the correlator (x 2 { x )x 2 (fy) when the couplings 52, ff3 and (74 
are neglected since the term gih\ 2 is a source for x 2 when the h kinetic term is omitted. 

After having taken proper care of the ultraviolet divergent terms the finite contributions from the diagrams in 
equation J7J) are 

2/ t n 2 rr gi T ( 0103 \ T ( g\\ 

m{T) = m - T l^ x ~^ m A 92 ~^)~l^ m { 29i -^) ' (9) 

which clearly shows that the terms in parenthesis (third and fourth terms) are subleading with respect to the second 
one near the transition point. Equation illustrates how the nearby phase transition is directly felt by the non-order 
parameter field, and furthermore, gives the general prediction that the screening mass of the singlet field must decrease 
close to the phase transition. If we stop the analysis here, at one loop level, we predict the following critical behavior 
for Am 2 — m 2 (T) — m 2 , where m is the mass at a temperature close to the critical point: 

Am 2 (T) = - JL*L~t-*, T<T C (10) 
107T m x 



Am 2 (r) = -jl^-rf, T>T C 



A 



where M x = \/2\m x \ cx |T — T c \ u / 2 . This one loop result breaks down at the transition point, due to the infrared 
singularity. However, since h is not the order parameter field, its correlation length (i.e. 1/m) is not expected to 
diverge at the phase transition. In the following we investigate the behavior of the h screening mass near the phase 
transition by going beyond the one-loop approximation. 



A. Healing the IR behavior 



When analyzing contributions beyond one loop order to the h two-point function, the number of diagrams, and 
distinct topologies one needs to consider proliferates. We select a subclass of diagrams, that heal the infrared diver- 
gences, while capturing the essential physical properties of the problem at hand. Besides, this subclass of diagrams 
we consider is known to be exact in the large N limit of a theory with O(N) symmetry 0, and the summation can be 
performed exactly. We note that a large N approximation might not be the best choice for investigating the details of 
the phase transition, however we show that it well reproduces the behavior of the screening mass of hadronic degrees 
of freedom near the deconfinement phase transition as given by lattice simulations. A finite result for the two-point 
function, both in the broken and unbroken phases will emerge, while the expression for the mass of h turns out to be 
continuous across T c . 

In the unbroken phase (i.e. T < T c ) it is known, that for a generic O(N) theory in the large N limit the following 
chain of bubble diagrams represents the leading contribution, 



-Q- +h OQ^+hooo^ + 



5 



These diagrams constitute a geometrical series. For the O(N) theory this resummation procedure is exact in the 
large N limit [j| . We choose the same set of diagrams in the unbroken phase. Denoting the dimensionless single loop 
integral by X, the expression for a single diagram in the chain made of p loops, and p — 1 vertices of x 4 -type is: 

|(_^)P-ilP. ( i 2) 
The exact sum of all the diagrams in the chain is simply: 



l 



q\X ^= , with 1 = —— . (13) 

1 1 + ^ 8ttto x 



Using eq. 1)13(1 the following expression for m follows: 



m 2 (T) = m 2 -r g ' (14) 

107T m x + A 1 

which is finite at T c , where it yields: 

m 2 (T c ) = m 2 -^. (15) 

This is the main result of Q, predicting that close to the phase transition the singlet state must have a decreasing 
mass parameter, associated with spatial correlations. More specifically, the drop at the phase transition point is 
given by the ratio of the square of the coupling constant governing the interaction of the singlet state with the order 
parameter (i.e. gi) to the order parameter field self-interaction coupling constant A. In this way, via the drop of the 
singlet field at the phase transition, one can derive a great deal of information about the phase transition, and about 
the order parameter itself. 

The analysis is much more complicated in the broken phase. Indeed, when T > T c \ develops an average value 
(x) — v i which induces one also for h as shown in eq.10). This is in agreement with the results found in 0]. So on 
general grounds, the fields x an d h mix in this phase. Defining the mass eigenstate fields as H and H, these are related 
to h and \ vi a: 



h \ _ f cos9 -sm6 \ ( H 
X J \ sin (9 cos 6 



(16) 



The mixing angle 9 is proportional to giv/m 2 , and therefore, the mixing can be neglected within the present 
approximations, simplifying the analysis considerably. Like for the symmetric phase, we consider only the effects due 
to the x loops for the h propagator. Due to symmetry breaking we now need to consider also the trilinear x coupling 

-^-X 3 , (17) 

which is expected to affect the analysis. At one loop level the only diagram to compute is again the one in eq. JSJl, 
with m x replaced by M x . Hence, on general grounds, we predict a drop in the mass of h also on the right hand side 
of T c . In this phase the infrared divergence is still present at one loop level. Curing such divergence is now more 
involved due to symmetry breaking. In the case of the large N limit of O(N) symmetry, one can show, that diagrams 
with trilinear vertices are again suppressed relative to the simple bubble diagrams. 

Here we go beyond the large N limit by computing a new set of diagrams, which can be evaluated exactly, and thus 
capture relevant corrections due to symmetry breaking, neglected in the large N limit. The new chain of diagrams 
we compute has terms of the form: 

-00 • • 00 • • <D- 

This class of diagrams has knowledge about the onset of the symmetry breaking via the presence of the trilinear 
vertices, and, in that respect, is the simplest extension of the chain of simple bubble diagrams. Another amusing 
property of (|18f) is, that the sum can be performed exactly, as we now demonstrate. 

Denoting the number of simple bubbles by p, and the number of bubbles with trilinear vertices by q, the total 
number of loops is L = p + 2q, and the number of x 4 vertices is p + q — 1, while that of the trilinear vertices is 2q. 
The loop integral of the single bubble is I. The loop integral of the bubble with two trilinear vertices yields X 2 /9. 
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The generic single multi-loop diagram of type with p ordinary bubbles and q bubbles with two trilinear vertices, 
yields the following contribution: 



4 V 4 v V 3! / \ Xv2 J 329 V 2 / \ 1 J 2 v 2 y V 3 / v 3 



(19) 



On the left hand side of the equation the binomial factor corresponds to the number of different ways to distribute q 
rungs into p + q loops, and the preceding factorials are the associated combinatorial factors, due to the permutation 
of the vertices. 

The infinite sums over q and p are evaluated by replacing the summation over p with the summation over the total 
number of loops L. For fixed L, q can take values q = 0, 1, . . . , [L/2\ , where [^J is the largest positive integer smaller 
than or equal to x. Then the sum becomes 

where the function Ul{x) appearing in the intermediate step is the Chebychev polynomial of the second kind. Re- 
markably, the final sum can be performed exactly, yielding for T > T c again an infrared finite result: 

m 2 (T) = m 2 - 4? - , (21) 



2 i + |l+4ix2 



T 



8ttM x 



(22) 



We see that m{T c ) from the broken side of the transition equals exactly the one from the unbroken side of the 
transition, even when departing from the large N limit. The mass squared of h is a continuous function through 
the phase transition, and the associated correlation length remains finite. This result does not hold order by order 
in the loop expansion, but only when the infinite sum of the diagrams is performed. Although admittedly, there is 
ambiguity in selecting the class of diagrams to take into account, guided by the large TV limit in the symmetric phase, 
we considered only the bubble chain of diagrams. However, in the broken phase we chose to resum a new and richer 
class of diagrams. The new class contains not only the relevant large N contribution, but also carries information 
related to the spontaneously broken symmetry. This class, although more involved, can be summed exactly. 

In order to disentangle relevant properties of the phase transition, we construct slope parameters for the singlet 
field: 

^+ 1 dm 2 (T) 

V ± = lim - — — , (23) 

with Am 2 (T c ) — g\j\ . The functional form of T> + and T>~ is the same, provided that the same class of diagrams is 
resummed on both sides of the transition. In this case we have: 

V- = lim — (24) 

XT dT ^ t - dT K ' 

in the symmetric phase, and to obtain T> + in the broken phase is sufficient to replace m x with M x — \ 2\m x \. While 
the mass of h remains finite at T c , its slope encodes the critical behavior of the theory. For example, if m 2 vanishes as 
(T c — T) v close to the phase transition (with the correlation length £ oc \T — T c |~"/ 2 ), then scales with exponent 
{vjl — 1). A difference in the functional form of the slopes may emerge, when on the two sides of the transition 
different class of diagrams are resummed. Using the wider class of diagrams in the broken phase considered above 
while retaining just the simple bubble sum in the unbroken we determine: 



V+c±-3^\m x \V-. (25) 



This is due to the onset of spontaneous symmetry breaking, i.e. in the broken phase the class of resummed diagrams 
contains trilinear type of interactions. We identify thus a less singular behavior with respect to the simple sum of 
bubbles. More specifically, the scaling exponent for T> + is now (v — 1). While the explicit relations between the 
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FIG. 2: Left panel: behavior of the mass of the singlet field close to the phase transition as function of the temperature. Right 
panel: Lattice data; figure from 

scaling exponents and the slopes are only valid within the given summation scheme, these quantities are, nevertheless, 
a measure of the critical behavior near the phase transition. Only experimental results will be able to select which 
class best describes the data. 

In figure|21we schematically represent the behavior of the h screening mass as a function of the temperature, in units 
of the critical temperature, for m 2 , cx (T c — T) . The left panel of fig. |21 illustrates the rapid decrease of the singlet 
field screening mass in the critical region. The right panel of fig. [21 is the lattice data of 6] obtained for the three 
color Yang-Mills theory. The mass in |(| is the glueball screening mass. The resemblance between our results and the 
lattice results is intriguing, but more lattice data is needed in order to distinguish between the above possibilities, 
and quantitatively determine the size of the drop. 

It is common practice to isolate the order parameter field, and more generally, the light degrees of freedom, since 
they are expected to be the relevant degrees at low energies. While this is certainly correct, we have explicitly shown 
how the essential features of the phase transition are encoded in the non-order parameter fields of the theory as well. 
Our analysis is useful whenever the order parameter field is neither a physical quantity, nor is phenomenologically 
accessible (e.g. the Polyakov Loop). 

IV. TIME DEPENDENT ORDER PARAMETER FIELD 

In this section we investigate the effects of the four dimensional x — x( x j t) order parameter field on the properties 
of a singlet h = /i(x, t) field. Some of these effects turn out to be significantly different than for the time independent 
case. The order parameter field is now a physical field, and as such, it can propagate in time, and can be canonically 
quantized. The four-dimensional Lagrangian of the renormalizable theory we use to define our Feynman rules is: 

U = -d.hd^h + -d^X - —h 2 ~ ^ X 2 " V ~ - h X 2 ~ 9 - h 2 X 2 - ^h 3 - ^ . (26) 
2 2 M 2 2 4! 2 4 3! 4! v ' 

All coupling constants are real, and A > 0, 54 > 0, and gi are dimensionless, while g\ and (73 have mass dimension 
one. The largest scale is the mass m of the h field, which, by our assumption, is again much larger than m x , the 
mass of the x field, and it is also much larger than the temperatures involved. We write then the couplings in terms 
of dimensionless ones via g\ — g\m , gi — §2 , and g% = g^m . The Lagrangian (|26[) contains all the relevant and 
marginal operators: gihx 2 and g^h 3 are relevant operators, while, for example, gi^X 2 is marginal. We expect the 
relevant physics to be well described by this set of operators. 

The field x is subject to Zi symmetry. Having the Yang-Mills deconfining phase transition in mind, in the theory 
with static order parameter field we chose the symmetric phase to be at low temperatures, and thus the symmetry 
was broken for temperatures greater than a phase transition temperature. Accordingly, m 2 , > for T < T c and 
m 2 , < for T > T c . As we anticipated though in the previous section, when the order parameter field is non-static we 
cannot simply choose the direction of symmetry restoration a priori. In four-dimensions a Bose-Einstein distribution 
emerges for x, and thus thermal fluctuations become important. Thermal fluctuations have the tendency to restore 
symmetry, destroying the possibility for the formation and existence of any physical condensate at high temperatures. 
The direction of the phase transition is therefore important, and is completely determined by the theory. In the 
following, we adopt the standard picture for the direction of the phase transition, consistent with the effect of higher 
order corrections to m x : symmetry restoration sets in at high temperatures p7j. This is one of the major differences 
between the theories with three- and four- dimensional order parameter field. 
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A. Diagramatics 



Based on Lagrangian (|26|l the diagrams contributing to the h two-point function at one-loop order are shown on 
(|7|). Even though these are the same as in the three dimensional theory, here they are computed taking into account 
the time dependence of \- We can show, that finite temperature corrections from all diagrams that involve an h loop 
are Boltzman suppressed due to the heaviness of h . Accordingly, the relevant diagrams in Q for the singlet h field 
are those on the first line: the first (tadpole A), the second (eye), and the third (tadpole B). Unlike in the three 
dimensional, time independent case, where UV divergent tadpoles have been removed by renormalization, in the four 
dimensional, time dependent case, these diagrams also provide non-negligible temperature contributions. Recall, that 
previously only the zero mode ft-o coupled in the /ix 2 -interaction. This is clearly not the situation here: all modes of h, 
and all Matsubara modes of x participate in the interactions. The diagrams are evaluated using standard techniques 
of the imaginary time formalism Q . 



d 3 k 1 52 f d 3 k 1 

(2tt) 3 uo 2 n + uj 2 ~ J (2tt) 3 277 

where w = \/k 2 + m x 2 is the energy, and uj m = 2nmT is the Matsubara frequency of the internal x-line. The second 
line in (|27|) is the result of the frequency summation, and f(x) — l/(e x — 1) is the Bose-Einstein distribution function. 
The term proportional to the number 1 before the distribution function represents the vacuum contribution, and 
is divergent in the ultraviolet. This divergence is absorbed in the mass renormalization process. The second term 
contains no divergence, since it is regularized by the distribution function. We focus on the physics in the phase 
transition region, where m x <C T. Since the mass of the order parameter field, m x , vanishes at T c , we use the high 
temperature limit to obtain an analytic result for the h self-energy contribution 

^-tadpole — ■ (^) 

Similarly, the other tadpole contributes with 

YiB 5i53 (9Q , 

^tadpole - 24 TO 2 ' 

These contributions are real and provide temperature dependence to the h mass, not present in the three dimensional 
theory. In terms of dimensionless couplings 

T 2 

^tadpole = (92 - gih) 24 ■ ( 30 ) 

Contribution to the self-energy of h from the eye diagram is 

d 3 k 1 1 



5i 



) 3 LU 2 n + LU 2 (LO m - LU n ) 2 + LO 2 

2 r d 3 k i i 



(2tt) 3 u 4lu 2 + uj 2 



(l + 2/( W /T)). (31) 



Here oj n is the Matsubara frequency of the external line and ui = y/k 2 + m x 2 . We have restricted the study to zero 
external momentum. In order to obtain physical quantities from this real Euclidean integral we analytically continue 
this to Minkowski space. After the replacement iuj n —> E + it we obtain 

IW*) = -2(|)7 (32) 

Another significant difference with respect to the theory with time independent order parameter field is the existence 
of a real and an imaginary part of the self-energy, that contributes to the two-point function \E 2 — m 2 — TL(E)\ 
Here m is the tree level mass of h and at one-loop order n = Tltadpoie + n e j, e . The real and imaginary parts can be 
extracted from by using 

' P-=fin6(z). (33) 



z ± ie z 
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Accordingly, the real part is given by the principal value of (13211 

Re U eye (E) = 2 (f ) 2 J^P / d|fc£_i_(i + 2f(u/T)) . (34) 

The term proportional to the number 1 before the distribution function is again the vacuum contribution, divergent 
in the ultraviolet, and is removed by usual renormalization. This real part represents a shift in the mass squared of 
the h, and for h at rest, through its definition, determines the pole mass M: 

M 2 - to 2 - U tadpole - Re U eye (E = M) = . (35) 

Self-consistent numerical solutions of the above equation show that the large tree level mass is dominant, and loop 
corrections are negligibly small in the temperature range of interest, i.e. near the phase transition. Another interesting 
fact is, that M ~ m acts as an infrared cutoff guaranteeing the absence of infrared divergence for the pole mass. As 
a consequence, in this case there is no need for the resummation of higher order diagrams. An analytic expression in 
powers of T 2 /m 2 , and setting m x = , can thus be obtained in the region we work in, namely m x <C T <C m, 

(9i\ 2 f^T 2 4vr 2 r 4 \ g\ 7 T 2 

Re IL eve (E = m = -2 ( — - + r ~ -^-m 2 — . (36) 

v K ' 12/ V 6m 2 15 m 4 7 12 m 2 y ' 

The absence of the infrared divergence in the pole mass is another relevant difference between the two theories 
discussed in this paper. This becomes apparent when summing all the one- loop contributions from l|3U|l and l|34|l : 



M = m 



1 + (.92 - 51.93 - 2g 2 ) -^-j + O ( 

24 ttt/ V ni^ 



(37) 



The imaginary part of the self energy represents the net rate of /i-decay, and its evaluation provides no difficulties: 



/oi\2 1 / 4m 2 

Im U eye (E) = -2 (!) —e(E-2m x )^l- —±(l + 2f(E/2T)). (38) 

By rewriting (1 + 2/) = (1 + f) 2 — f 2 it is easy to understand that above the threshold, E = 2m x , (|38|l determines 
the decay process h — > XX > an( i the annihilation \X ~* h , weighted by the thermal distribution, with the net decay 
rate given by 

„ Im n p „p (E) 

r = ■ (39) 

Unlike in the three dimensional case, here there are finite temperature corrections also to the two-point function of 
the now dynamic order parameter field. The full set of one-loop diagrams contributing to the x field in the symmetric 
phase, together with their corresponding combinatorial factors is: 



=0= = 12=0= + 2=o= + 4 q + 2 9 +« 9 

(40) 

Due to the existence of a cubic x self-coupling in the broken phase, there are, in addition to 1)41)0. two more diagrams: 

(41) 

Here we are not going to discuss further details about the x two-point function. The evaluation of these diagrams 
is standard. There are two points we wish to stress though: First, that contrary to the three dimensional theory, 
there are explicit thermal fluctuations providing temperature dependence to the mass of the order parameter field; 
Second, this mass is zero at the critical temperature of the second order phase transition, and has thus m 2 oc (T c — T) 
behavior. 
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B. Static Properties 

In the three-dimensional theory we found an infrared divergent result for the one-loop contribution to the h two- 
point function. This divergence came from what we call the eye diagram. In four dimensions, as we have shown above, 
besides the eye diagram there are two tadpoles also contributing. Here we show that the static limit of the two-point 
h function for the time dependent order parameter field displays the same features as for the time independent order 
parameter field case. 

Here it is important to distinguish between the pole mass and the screening mass. The pole mass, M, is defined 
through (|35|l as the pole in the full two-point function. The screening mass, m s , is defined by the location of the pole 
in the static propagator for complex momentum p = im s , 



P 



U{E = 0,p) = 0. 



(42) 



In the small momentum limit this leads to the following definition 



lim U(E = Q,p) 
p— >o 



(43) 



In three dimensions the pole and screening masses are one and the same thing, since the static propagator is the full 
propagator. 

We have shown, that when looking at the pole mass the IR problem of the eye diagram is regulated by the heavy 
h mass. When analyzing the screening mass, however, the one-loop IR divergence present in the three dimensional 
case is recovered. Based on the above definition, we set E = in expression i|34|) . This corresponds to effectively 
reducing the 4 dimensional theory to a 3 dimensional one. To display the relevant infrared contribution we take the 
high temperature expansion in l(34j) which yields: 



9i\ 2 T_ 
At: 2 



dp 



(P 2 



,2\2 



T 



327T771 



(44) 



x 



Hence, for the static limit the phase transition region is dominated by the same type of infrared divergence we 
encountered in the time independent order parameter field case. Note, however, that the numerical constant in eq. 
(14411 differs from that in the second term of eq. © ■ The reason for this is that the reduction, both here and in section 
IIIII was done only for the modes of h, and thus in eq. 144(1 . contrary to 10, all Matsubara modes of x contribute. 
By combining all the diagrams we find for the screening mass at one-loop order 



1 



9l 



T 



327T m y 



{fa fa - fa) 



T 2 



24 m 2 



(45) 



showing clearly the eye contribution as the infrared dominant one. Note, that above we have tree level coefficients. A 
complete investigation would require renormalization group analysis, but this is beyond the scope of this paper. We 
thus recovered our universal behavior from before. 



V. CONCLUSIONS 



When analyzing many physical situations, it is common practice to isolate the order parameter field, and more 
generally, the light degrees of freedom, since these are expected to be the relevant states at low energies. While this 
procedure certainly is reasonable, in nature most of the physical fields are neither order parameter fields, nor light at 
all. In order to extract information from these heavy states, we needed first to determine new and universal features 
associated with them. We have used a general strategy proposed first in 0, Q , according to which we couple the light 
degrees/order parameter fields to the heavy fields in the most general way, and then truncate the theory by retaining 
all the relevant and marginal operators in the Lagrangian. In doing so, the theory is fully renormalizable while 
capturing the relevant contributions. The operator set is further constrained by imposing all the relevant symmetries 
of the problem at hand. In order for our procedure to work, we also assume, that the other physical states of the 
theory have masses larger than our non-order parameter field. In this way we can, formally, integrate these states 
out, and their effects are absorbed in the modified couplings of our effective Lagrangian. 

For both the time independent and time dependent order parameter field we have shown, that the spatial correlators 
of the non-order parameter field are infrared dominated, and hence can be used to determine the properties of 
the phase transition. We have determined the general behavior of the screening mass of a generic singlet field, 
and have shown how to extract all the relevant information from such a quantity. Somewhat surprisingly, we have 
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demonstrated, that the pole mass of any non-order parameter physical field is not infrared dominated. Our results can 
be immediately applied to any generic phase transition. We have used as relevant example, for the time independent 
order parameter field case, the deconfining transition of Yang-Mills theories as also explained in 0, Lattice 
simulations Q for the deconfining phase transition support our results. Different a ppr oaches have been used in 
literature to study the deconfining phase transition 0, [iS [13, [13, EE Q O EE [O, 111, E3, EL . Here we 
have been able to unify some of them. In fact, by demonstrating that there exists an extended universal behavior 
for singlet fields, we have uncovered the relation between phenomenologically oriented models which use the glueball 
Lagrangian [TH Li EE EE [HI to describe the deconfining phase transition, and the ones using the symmetries of the 
Polyakov loops [H Ej, El H EE . 

The induced critical behavior is universal, however the quantitative details depend on the strength of the couplings 
between the fields as well as the resummation procedure. We have considered only one non-order parameter field, 
but many are expected to display a similar behavior. For the Yang-Mills deconfining phase transition lattice QCD 
simulations are able to determine the coupling strength of any glueball state to the Polyakov loop by following the 
temperature dependence of screening masses of such states. 

Our analysis suggests, that monitoring a number of spatial correlators, or more specifically their derivatives, is an 
efficient and sufficient way to experimentally uncover the chiral/deconfining phase transition and its features. 
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